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Abstract: We investigate the deconfining phase transition in SU(3) pure gauge theory 
and in full QCD with two flavors of staggered fermions. The phase transition is detected 
by measuring the free energy in presence of an abelian monopole background field. In the 
pure gauge case our finite size scaling analysis is in agreement with the well known presence 
of a weak first order phase transition. In the case of 2 flavors full QCD we find, using the 
standard pure gauge and staggered fermion actions, that the phase transition is consistent 
with weak first order, contrary to the expectation of a crossover for not too large quark 
masses and in agreement with results obtained by the Pisa group. 
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1. Introduction 

Understanding QCD thermodynamics is one of the most intriguing issues of contemporary 
physics. Indeed the study of QCD at high temperature (and density) [1-4] is relevant for 
both high energy physics (e.g. ultrarelativistic heavy ion collisions) and astrophysics (com- 
pact stars). Moreover addressing QCD thermodynamics could shed light on the problem 
of color confinement and chiral symmetry breaking. 

To detect the deconfinement phase transition in pure gauge theories the expectation 
value of the trace of the Polyakov loop is commonly used as an order parameter. In 
presence of dynamical fermions the Polyakov loop ceases to be an order parameter since 
Z(N) symmetry is no longer a symmetry of the action. Alternatively, in order to determine 
the finite temperature phase transition, one usually studies the chiral condensate, which 
however is not related to confinement, but to chiral symmetry breaking, and in its turn is 
not a good order parameter at non-zero quark masses. 

A mechanism for color confinement based on dual superconductivity of the QCD vac- 
uum by abelian monopole condensation has been proposed in [5-7]. A disorder parameter 
which is related to abelian monopole condensation in the dual superconductivity picture of 
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confinement has been developed by the Pisa group and consists in the vacuum expectation 
value of a magnetically charged operator, (fi). In Refs. [8-10] it has been shown that (/x) is 
different from zero in the confined phase of SU(2) and SU(3) pure gauge theories, that it 
goes to zero at the deconfining phase transition, and that this is independent of the abelian 
projection chosen to define the magnetic charge. The same has also been verified in the 
case of full QCD [11]. 

In the present paper abelian monopole condensation is detected by looking at the 
free energy [12, 13] in presence of an abelian monopole background field, which in turn 
is evaluated in terms of a gauge invariant lattice effective action. This lattice effective 
action is defined at zero temperature by means of the lattice Schrodinger functional and 
at finite temperature by means of a thermal partition functional, employed so far to study 
the vacuum dynamics of SU(3) lattice gauge theory at finite temperature in presence of a 
constant abelian background field [14] or in presence of an abelian monopole background 
field [12,13,15]. 

Since the free energy is related to the vacuum dynamics and not, like the trace of the 
Polyakov loop, to a symmetry of the gauge action which is washed out by the presence of 
dynamical fermions, we feel that it could be used to detect the finite temperature phase 
transition also in full QCD, as well as the order parameter (/j) employed in Ref. [11]. 

In the present paper we present a finite size scaling study of the free energy in presence 
of an abelian monopole background field both for SU(3) pure gauge and two-flavor QCD. 



The plan of the paper is the following. In Sections LI and L2 we introduce the 
lattice effective action at zero temperature and at finite temperature. In Section |1.3| we 
modify the thermal partition functional for including dynamical fermions. In Section || we 
apply our method to the case of SU(3) pure gauge theory and, by means of a finite size 
scaling analysis, we find that (as expected) the deconfining phase transition is first order. 
In Section ^ we explore QCD thermodynamics with two staggered dynamical fermions of 
equal masses using the standard staggered fermion action. A finite size scaling analysis 
suggests that the deconfining phase transition is "weak first order" , in agreement with the 
indications presented in [11, 16, 17]. In Section |] we sketch our conclusions. Finally in 
Appendix [A| we report our finite size scaling investigation of the U(l) lattice gauge theory. 
A preliminary account of the results of the present paper appeared in Refs. [18, 19]. 

1.1 The lattice effective action: T = 

In order to investigate vacuum structure of lattice gauge theories at zero temperature a 
lattice gauge invariant effective action r[j4 ext ] for an external background field A ext was 
introduced in Refs. [20,21]. It is defined as 

where Lt is the lattice size in time direction. A ext (x) is the continuum gauge potential for 
the external static background field, the corresponding lattice links are 

U% xt (x) = Pexp I dt iagAf\x + tak) (1.2) 
Jo 
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(a: lattice spacing, g: bare gauge coupling, P: path ordering operator). Z[A ext ] is the 
lattice partition functional 

Z[A cxt ] = [ VU e~ Sw , (1.3) 

JU k {S,xt=Q)=U^{g) 

Sw is the Wilson action and the functional integration is performed over the lattice links, 
but constraining the spatial links belonging to a given time slice (say xt = 0) to be 

U k (x,x t = 0) = U^(x) , (k = 1,2,3) , (1.4) 



U% xt (x) being the lattice version (see Eq. ( |1.2| )) of the external continuum gauge field 
A ext (x) = A ext (x)X a /2 (with A a /2 Gell-Mann matrices). The temporal links are not con- 
strained. In the case of a static background field which does not vanish at infinity we must 
also impose that spatial links exiting from sites belonging to the spatial boundaries (for 
each time slice Xt ^ 0) are fixed according to Eq. ( |1.4j ): in the continuum this last condition 
amounts to the requirement that fluctuations over the background field vanish at infinity. 

The partition function defined in Eq. (|1.3| ) is also known as lattice Schrodinger func- 
tional [22,23] and in the continuum corresponds to the Feynman kernel [24,25]. Note that, 
at variance with the usual formulation of the lattice Schrodinger functional [22,23] where a 
lattice cylindrical geometry is adopted, our lattice has an hypertoroidal geometry so that 
Sw in Eq. (1.3) is allowed to be the standard Wilson action. 



The lattice effective action T[^4 ext ] corresponds to the vacuum energy, So[j4 ext ], in 
presence of the background field with respect to the vacuum energy, £b[0], with A ext = 

T[A cxt ] — ► E [A ext ] - E [0] . (1.5) 

The relation above is true by letting the temporal lattice size Lt — > oo, on finite lattices 
this amounts to have Lt sufficiently large to single out the ground state contribution to the 
energy. 

Since the lattice effective action Eq. ( |l.ip is given in terms of the lattice Schrodinger 
functional, which is invariant for time-independent gauge transformation of the background 
field [22,23], it is as well gauge invariant. 

1.2 The thermal partition functional 

If we now consider the gauge theory at finite temperature T = l/(aLt) in presence of an 
external background field, the relevant quantity turns out to be the free energy functional 
defined as 

1 J L t \ Zt[0] J 



1.6) 



Zq-[A ext ] is the thermal partition functional [26] in presence of the background field A ext , 
and is defined as 



A 



cxt 



VUe~ Sw , (1.7) 

U k (x,L t )=U k {Sfi)=Ul^{x) 
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In Eq. (|1,7| ), as in Eq. (jOj), the spatial links belonging to the time slice xt = are 
constrained to the value of the external background field, the temporal links are not con- 
strained. On a lattice with finite spatial extension we also usually impose that the links at 
the spatial boundaries are fixed according to boundary conditions Eq. (1.4), apart from the 



case in which the external background field vanishes at spatial infinity (as happens for the 
monopole field), where the choice of periodic boundary conditions in the spatial direction 
is equivalent to Eq. (1.4) in the thermodynamical limit. If the physical temperature is 



sent to zero, the thermal functional Eq. ( |1.7[) reduces to the zero-temperature Schrodinger 
functional Eq. (fy|). 

The free energy functional Eq. ( |1.6[ ) corresponds to the free energy, i ? [yl ext ], in presence 
of the external background field evaluated with respect to the free energy, F[0], with 
A ext = 0. When the physical temperature is sent to zero the free energy functional reduces 
to the vacuum energy functional Eq. ( |1,1|) . 

1.3 Including fermions 

When including dynamical fermions, the thermal partition functional in presence of a static 
external background gauge field, Eq. (|1.7|) , becomes: 



A 



cxt 



VUVipVi>e-( Sw+SF) 

U k (L t ,x)=U k (0,x)=U c k xt (x) 

VU e~ Sw det M , (1.8) 

U k (L t ,x)=U k (0,x)=U^(x) 

where Sw is the Wilson action, Sf is the fermionic action and M is the fermionic matrix. 
Notice that the fermionic fields are not constrained and the integration constraint is only 
relative to the gauge fields: this leads, as in the usual QCD partition function, to the 
appearance of the gauge invariant fermionic determinant after integration on the fermionic 
fields. As usual we impose on fermionic fields periodic boundary conditions in the spatial 
directions and antiperiodic boundary conditions in the temporal direction. 

1.4 The monopole free energy and the monopole condensation 



We use the free energy functional Eq. Ql.6| ) to evaluate the free energy in presence of an 
abelian monopole background field (see Sect. |2.1| ). If there is monopole condensation, the 
free energy with the monopole background field is the same as the free energy without the 
monopole background field (^4 ext = 0). In other words it does not cost energy to create a 
monopole. Therefore by evaluating the free energy by means of Eq. ( |l.6|) we are able to 
detect abelian monopole condensation. 



Note that, as discussed in Sect. 1.1, our free energy functional is gauge invariant for 
time-independent gauge transformations of the external background field. This implies 
that we do not need to do any gauge fixing to perform the abelian projection. Indeed, after 
choosing the abelian direction, needed to define the abelian monopole field through the 
abelian projection, due to gauge invariance of Schrodinger functional for transformations 
of background field, our results do not depend on the selected abelian direction, which, 
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actually, can be varied by a gauge transformation. In this sense our definition is analogous 
to the definition given in Ref. [10], where the monopole field is defined without gauge fixing. 

To evaluate the monopole free energy by means of Eq. (|0|) we should compute the 
ratio of two partition functions. In order to avoid the problem of dealing with partition 
functions, we will compute instead [27,28] the derivative of the free energy F'(f3) = dJ-/d(3 
with respect to the gauge coupling (3 ((3 = 6/g 2 ). It is easy to see that F'((3) is given by 
the difference between the average plaquette < U^ u > obtained in turn from configurations 
with A ext = and A cxt / 



F'((3) 



d/3 



V 



(1-9) 



where V is the spatial volume. In Eq. ( |1.9| ) the dependence of the free energy functional 
J- on f3 (at fixed external gauge potential ^4 ext ) has been made explicit. Eventually, since 
J- = at (3 = 0, the free energy can be obtained by numerical integration of F'(fl) 



HP) 



F'(f3') d(3' . 



(1.10) 



2. SU(3) pure gauge 

2.1 The monopole background field 

For SU(3) gauge theory the maximal abelian group is U(l)xU(l), therefore we may in- 
troduce two independent types of abelian monopoles using respectively the Gell-Mann 
matrices A3 and Ag or their linear combinations. 

In the following we shall consider the abelian monopole field related to the A3 diagonal 
generator. In the continuum the abelian monopole field is given by 



gb a (x) 



a 3 "mon 



x x n 



\x\(\x\ 



x ■ n) 



(2.1) 



where n is the direction of the Dirac string and, according to the Dirac quantization 
condition, n mon is an integer. The lattice links corresponding to the abelian monopole field 
Eq. (|2.1| ) are (we choose n = £3) 



with 6™ Il (x) defined as 



1,2 

3 (2) 



e I. 2 



m 

e - i6 ^W 
1 



(2.2) 



(x 2 ~ X 2 ) 



1 



+- 



l^monl l^n 

(xi - X{) 

^-mori \Xr\ 



(x 3 - X 3 ) 
1 

(X 3 " X 3 ) 



(2.3) 
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Figure 1: F' vs. [3 (n m < 
(for the significance of spatial "fixed boundary conditions" see Section 2.1) 



10) for L s = 16, 24, 32 and Lt = 4. Spatial "fixed boundary conditions" 



where (Xi, X2, X3) are the monopole coordinates, x mon = (x — X). 

The monopole background field is introduced by constraining (see Eq. (|2.2| )) the spatial 
links exiting from the sites at the boundary of the time slice xt = 0. For what concern 
spatial links exiting from sites at the boundary of other time slices (xt 7^ 0) we consider 
two possibilities. In the first one we constrain these links according to Eq. ( |2.2| ) (in the 
following we refer to this possibility as "fixed boundary conditions"). In the second one we 
do not impose the constraint Eq. (2.2) on the above mentioned links (this possibility will 
be referred as "periodic boundary conditions"). 

We simulate pure SU(3) lattice gauge theory with Wilson action. The lattice geometry 
is hypertoroidal. The path integral is constrained as given in Eq. ( |2.2| ). It is important 
to observe here that since the background field potential is not an integration variable 
(i.e. it does not correspond to a thermalized field) it is not necessarily subject to periodic 
boundary conditions. 

The simulations were performed on lattices of different spatial sizes taking fixed the 
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Figure 2: F' vs. f3 (n mon = 10) for L s = 16, 24, 32 and L t — 4. Spatial "periodic boundary 



conditions" (for the significance of spatial "periodic boundary conditions" see Section 2.1). 



temporal extent (Lt = 4). We consider lattices with spatial volumes 16 3 , 24 3 , and 32 3 . 
Simulations were performed in part on a APE100/QH1 and in part on APEmille/crate in 
Bari. To upgrade SU(3) matrices we alternate a Cabibbo-Marinari heat bath sweep with 
one (or more) overrelaxed sweep. The statistics collected after 2,000 thermalization sweeps 
amounts at between 5,000 and 10,000 configurations for each value of the gauge coupling 
(5. The statistical analysis has been done by means of jackknife resampling. 

In Fig. |l]we report our numerical results for F' (Eq. (|1.S|)) versus for three different 
lattice sizes, n mon = 10, and spatial "fixed boundary conditions". The data display a sharp 
peak that increases by increasing the lattice spatial volume. The data reported in Fig. || 
refer to spatial "periodic boundary conditions" and display the same qualitative behavior, 
though with an increased signal in the peak region. This is to be expected, since the 
effective volume is larger than in the previous case with fixed spatial boundary conditions. 

By inspecting Figure || it is evident that F'(f3) = in a finite range of /3, starting from 
(3 = and below the critical coupling, signaled by a peak in F'(0). Therefore F((3) = in 
a finite range below the critical coupling, above which the gauge system gets deconfined. 
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The vanishing of the free energy implies abelian monopole condensation. 

On the other hand, F'{f3) becomes different from zero and increases with the lattice 
spatial volume near the critical coupling, as expected in presence of a phase transition. 
Moreover in the weak coupling regime F'(J3) stays constant and almost independent from 
the spatial lattice volume 1 . This corresponds to the classical monopole energy which de- 
pends linearly on (3. It is clear that in the deconfined phase it costs a finite amount of 
energy to create a monopole, and so, there is not abelian monopole condensation. 

In the following we will analyze the peak structure of F'{(3) near the critical coupling 
and show that the data display the right finite size scaling behavior with scaling parameters 
compatible with a first order phase transition. 



2.2 Finite Size Scaling 

As a first step we determine the value of the critical coupling (3 c {L e ^) in correspondence of 
each set of data (at different lattice sizes and different spatial boundary conditions). 
We find that our data for F'((3,L e ^) can be fitted according to 

n( T cS ) 

FULf)= W-l(LfT ' p<MLfl <2A) 

The quantity is defined as 

„ ff I L s — 2 for spatial "fixed b.c.'s" , 
Lf=\ 1 ' (2.5) 

[L s for spatial "periodic b.c.'s" . 

The definition of L| given in the previous Eq. ( |2.5| ) takes into account the circumstance 
that for spatial "fixed boundary conditions" the effective spatial volumes is reduced since 
the links exiting from the sites at the spatial boundaries of each time slice (xt ^ 0) are 
constrained to the value given in Eq. fl2.2| ). 

Unlike the case of U(l) lattice gauge theory, discussed in Appendix [A|, where a ~ 1, 
we find that our data are best fitted with a ~ 0.35. 

In Table [j] we display (5 c {L e ^) for both spatial periodic and fixed boundary condition. 
From the various results for f3 c {L e ^) we are able to get an estimate of /3 C = /3 c (L| ff = oo). 
We try the following fit 

P c (Lf)=p c + d 1 (Lf)- 1 ^. (2.6) 
We find a rather good fit (x 2 /d.o.f. ~ 0.8) with d± ~ 0.6 and 

C = 5.3228 ± 0.0024 , v = 0.334 ± 0.021 . (2.7) 

The critical exponent v is consistent with a first order phase transition where 1/v = 3. 



x We note a difference with the Pisa order parameter fi: in that case 4j ln(//) increases linearly with the 
spatial volume in the weak coupling limit, since [i is a magnetically charged operator subject to magnetic 
charge superselection in the deconfined phase [29]. 
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spatial "fixed boundary conditions" 


( spp Section \2 1 \ 




L cflf 


dr 


16 


14 


5 4033 (959) 


24 


22 


5 2976 f247l 


32 


30 


5.3138 (108) 




spatial "periodic boundary conditions 


" (see Section ^4|) 


L. 


Lf 




16 


16 


5.3296 (33) 


24 


24 


5.3200 (56) 


32 


32 


5.3199 (100) 



Table 1: The value of j3 c in correspondence of different lattice sizes and different spatial boundary 
conditions, obtained by fitting Eq. (2.4) to the lattice data for F'(j3, L e s s ) (see Eq. (1.9)). 



The behavior of F'(P, Lf) in the critical region can be investigated by using finite size 
scaling techniques. Equations (|2.4j) and (p.6Q suggest the following scaling law 



F'((3,L 



cff\ 



a i 



{Lf) a/v 



\P-Pc- di(Lf )-V»\ a \(Lfy/-(P - p c ) - di 
so that F'(/3, Lf)/(Lf) a l v is a universal function of the scaling variable 

x = {Lf) l l»(f3-f3 c ). 
Note that Eq. Q2.8]) gives a sensible result in the thermodynamical limit L 



(2. 



(2. 



eff 
s 



OC 



F'W) 



ai 



\P~Pc 



P<Pc 



(2.10) 



while F'(P = 0c, Lf) diverges like (Lf ) a l u when Lf oo. 
Accordingly we fitted our lattice data with the scaling law 



F'(P,L 



f) 



ai (Lfr 



\(Lfy/»{p-p c )-d x 



(2.11) 



where we expect that 7 = a/ v. 

The output of the fits are reported in Table |2[ The parameter d\ has been fixed at 
the value obtained with the fit Eq. ( |2.6| ). We see clearly that ajv agrees with 7 within 
statistical errors, as expected if we want a sensible result in the thermodynamical limit 
Lf 00 (see Eq. pUD). 

In Figs. |, I we plot F'(p,Lf) rescaled with the factor l/(Lf)^ versus the scaling 
variable x defined in Eq. ( |2.9| ), respectively for spatial "fixed boundary conditions" and for 
spatial "periodic boundary conditions" . The full line is 



F'(P,L 
(Lf). 



eff\ 



\{Lf) 1/v (P - Pc) - dA 



(2.12) 
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Figure 3: The values of (3 C (L C S S ) (open circles), versus 1/L C S S , obtained at different values of L e s s 
and different spatial boundary conditions (see Table [l]). L c s s is denned in Eq. (2.5). The dotted line 
is the fit Eq. (|j). The value /3 c {Lf = oo) is also plotted (full circle). 



We find that the scaling relation holds quite well for a very large range of x. The quality 
of the scaling can be inferred looking at Figs. ^| and ||. Moreover looking at Table ^ we 
remarkably see that the critical parameters a, /3 C , is, and 7 agree for both sets of data. 

Finally in Fig. ^ we display F'(/3, L| ff )/(L| ff ) 7 versus the scaling variable x for both 
periodic and fixed spatial boundary conditions. We see that the numerical data nicely 
display the same scaling behavior. 

It is interesting to comment on the behavior of exp(— F({3)/T), which is the analogous 
of disorder parameter of Refs. [8-10], in the thermodynamical limit implied by Eq. ( 2.10) ). 
Indeed we have that 



exp ("^) = exp (-±J F'((3') d^j , ft < P < 0c , (2.13) 
while we already know that exp(— F(0)/T) = lfor (3 < (3q irrespective of the lattice size 
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spatial "fixed boundary conditions" 


(see Section [2.1) 




I He v 


d\ OL 


12.199 


1.247 5.3251 0.335 


0.6 0.351 


±3.9004 


±0.089 ±0.0110 ±0.026 


constant ±0.035 




spatial "periodic boundary conditions 


" (see Section ) 


Oi 


7 0c v 


d\ a 


13.461 


1.510 5.3222 0.340 


0.6 0.347 


±1.337 


±0.555 ±0.0013 ±0.020 


constant ±0.009 



Table 2: The values of the parameters obtained by fitting Eq. (2.11) to the data for the derivative 



of the monopole free energy Eq. (1.9) on lattices with spatial volumes 16 3 , 24 3 , and 32 3 and spatial 
"fixed" or "periodic" boundary conditions respectively. 



(see Fig. (§)). From Eq. (|2.1C| ) we get 



exp 



F(J3) 



exp 



1 Ol 



T 1-a 



(\0o ~ 0c 



1 1-a 



\0-0cV- a ) 



0o<0<0c (2.14) 



So that F(0) decreases when — > /3 C if < a < 1 tending to a finite value at = C . 
On the other hand, for a = 1 it is easy to see that exp(— F(0)/T) decreases to zero as 
a power of (0 C — 0). Thus we conclude that for a < 1 we have a discontinuous jump of 
exp(— F(0)/T) at C and the strength of the discontinuity weakens when a — * 1. However, 
it must be stressed that the discontinuous jump of exp(— F(0)/T) is exceedingly small so 
that exp(— F([3)/T) decreases almost continuously toward zero when — ► /3 C . 

2.3 The deconfinement temperature 

Before ending the discussion of pure SU(3) gauge theory, we would like to remark that 
using the data for F'(0) we are able to get an estimate of the continuum extrapolated 
critical temperature T c in units of the square root of the string tension \fa. We will show 
here that our estimate is consistent, even though with a quite large error, with the updated 
value in the literature. 

In fact, by fitting Eq. (2.4) to our data for F'(0) obtained on lattices 32 3 x Lt (Lt = 



6,7,8, with spatial periodic boundary conditions), we are able to get an estimate of the 
critical coupling c {Lt) corresponding to each Lt. To fix a physical scale we consider the 
string tension a^/cr at each value of c {Lt). The string tension is obtained on a symmetric 
lattice with the Wilson action. To this purpose we use the data for the string tension as 
parameterized in Eq.(4.4) of Ref. [30]. 

In Fig. [7| our data for T c /\/a are displayed versus aT c . The continuum extrapolation 
using an ansatz quadratic in aT c gives the estimate T c /y/a = 0.635±0.147 to be confronted 
with the average estimate in the literature T c j ' \J~o = 0.640 ±0.015 (see Table 3 of Ref. [31]). 

Our conclusion is that our method allows to get an estimate of T c j \fa although the 
statistical uncertainty is quite large, mainly due to a large error in the evaluation of c (Lt = 
6). A better estimate of C and then of T c /y/a could be achieved by means of the density 
spectral method. 
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Figure 4: The derivative of the monopole free energy with respect to the gauge coupling f3 
(Eq. (|L9|)) rescaled with (Lf)^ versus the scaling variable x (Eq. ( |2.9| )). Data refer to simula- 
tions with spatial "hxed boundary conditions". The full line is the scaling curve Eq. (2.12) with 
the parameters from Table @. 



3. QCD with two dynamical flavors 



We will account for results achieved from the study of QCD with two dynamical fermions 
using the standard staggered fermion action. As discussed in the Introduction, we simulate 
the theory with a "cold" time-slice (say xt = 0) where the spatial links are constrained 



to be a (lattice) abelian monopole background field (Eq. (p.2[)). According to Section 2.1 



the spatial links exiting from the sites belonging to the boundaries of the Xt ^ temporal 
slices are also constrained according to Eq. (|2.2|): we have not considered the possibility of 
spatial periodic boundary conditions in this case. 

We compute the derivative of the monopole background field free energy with respect 
to the gauge coupling, as given in Eq. fll.SD , where now the expectation value is evaluated 



with the full QCD action (details and numerical results in Sections 3T and |3.2| ), Our main 
goal is to try to use our data for F'{(3) on different spatial volumes and different bare quark 
masses to infer the critical behavior of two flavors full QCD near the deconfining transition. 
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Figure 5: The derivative of the monopole free energy with respect to the gauge coupling f3 
(Eq. (|T|)) rescaled with (Lf)^ versus the scaling variable x (Eq. ( |2.9| )). Data refer to simula- 



tions with spatial "periodic boundary conditions". The full line is the scaling curve Eq. (2.12) with 
the parameters from Table @. 



3.1 Numerical simulations 

We used a slight modification of the standard HMC R-algorithm [32] for two degenerate 
flavors of staggered fermions with quark mass m q : the links which are frozen are not evolved 
during the molecular dynamics trajectory and the corresponding conjugate momenta are 
set to zero. We have collected about 2000 thermalized trajectories for each value of (3 at 
L s = 16, 20 and about 1000 thermalized trajectories for each value of (5 at L s = 32. Each 
trajectory consists of 125 molecular dynamics steps and has total length 1. The computer 
simulations have been performed on the APEmille crate. 



3.2 Numerical results 

In Fig. H we compare F'(/3) for two staggered degenerate flavors with quark mass m q = 
0.075 and in the quenched case. In the 2 flavors full QCD case the signal in the peak region 
gets enhanced with respect to the quenched case and the position of the peak shifts to a 
smaller value of (3. 
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Figure 6: We plot together data from Fig. [| (fbcs) and Fig. || (pbcs). 



In Fig. |we display F'{(3) in the case of 2 dynamical flavors (m q = 0.075, 16 3 x4 lattice) 
compared with the chiral condensate < ipip >. Data displayed in Fig. || suggest that the 
peak in F'{(3) corresponds to the drop of the chiral condensate. We varied the lattice size 
L s (L s = 16,20,32) and the staggered quark mass m q (m q = 0.075,0.2676,0.5003). At 
fixed m q , as in the quenched case, we find that the sharp peak increases with the lattice 
spatial volume. Moreover at fixed spatial volume the critical coupling depends on m q . 



3.3 Finite Size Scaling 

We would like to use the numerical data collected at different lattice size and staggered 
quark mass to perform a finite size scaling analysis. According to Eq. ( |2.8D we try the 
scaling law 

F'^Lf) = l/T ^" l{LfV r r ria , (3.1) 



\(Lf)^(j3-p e (m q ))-d 



where the critical coupling (3 c {m q ) depends on the quark mass m q . The dependence of the 
critical coupling (3 c {m q ) on the quark mass is determined by the chiral critical point [33-35]. 
In the thermodynamical limit, by known universality arguments the critical couplings will 
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Figure 7: Continuum extrapolation for T c /*Ja (full circles). The extrapolated value is T c j^fa 
0.635 ±0.147 (full square). 



scale like 



f3 c (m q ) = (3 c (m q = 0) + cm, 



1/(35 



(3.2) 



where 1//35 is a combination of critical exponents which for the case of 2-flavors QCD are 
expected to be those of the three-dimensional 0(4) symmetric spin models: 



-0.5415. 

pd 

Inserting Eq. ( |3.2| ) into Eq. (|3.1| ) we are lead to the following scaling law 



F'(P,Lf,m q ) 



\(Lf)V"(p-p c (p)-cn%)-d 1 \ 



(3.3) 



(3.4) 



where again 7 = a/v assures a sensible thermodynamical limit. Note that, to take care of 
finite volume effects, the exponent rj is expected to be: 



v pd 



(3.5) 
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Figure 8: F' vs. (3 in the case of pure SU(3) gauge theory (open circles) and full QCD with 2 
dynamical flavors (open squares). Lattice size is 32 3 x 4. 



where u', (3, and 5 are the chiral critical exponents. 

Indeed, Eqs. flOp and (|3.5p assure that in the scaling region 



In our case the relevant chiral critical exponents are those of the three-dimensional 0(4) 
symmetric spin models where [34] 

z/ = 0.7423, v c = 0.4019. (3.7) 



In Table y we report the results obtained by fitting Eq. ( |3.4| ) to all our lattice data. 
From Table |3| we can see that a/v = 2.35 ± 0.34 consistent with 7 = 2.00 ± 0.47 (see 
Section |2.2[) . Concerning the parameter rj we find that it is poorly determined by our data. 
If we constrain rj in our fit we get 77 = 1.10 ± 0.19 which, together with v = 0.31 ± 0.03 
leads to (see Eq. (^5|)) v c = 0.34 ± 0.07 consistent with the value reported in Eq. (|3~^). 
However if we release the constraint on 77 our data can also be fitted with smaller values 
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Figure 9: F' on a 16 3 x 4 lattice in the case of full QCD with 2 dynamical flavors, is displayed 
together with chiral condensate. 





spatial "fixed boundary conditions" 


(see Section 2T|) 




dl 


7 &(0) c rj 


v d\ 


a 


79.4 


2.00 4.9933 0.54 1.10 


0.31 0.6 


0.728 


±76.6 


±0.47 ±0.0138 ±0.11 ±0.19 


±0.03 constant 


±0.078 



Table 3: The values of the parameters obtained by fitting Eq. (3^4) to the data for the derivative 
of the monopole free energy Eq. (L9) in two-flavors full QCD on lattices with spatial volumes 16 3 , 
24 3 , and 32 3 and L t = 4. 



for rj without altering significantly the other parameters. Moreover by confronting the 
exponent a in Table |3] with the corresponding value for the SU(3) pure gauge in Table || 
we conclude that our data for full QCD with two dynamical flavors are compatible with 
a first order phase transition {y = 0.31 ± 0.03) but, in the sense discussed at the end of 
Section 2^2, this is weaker than in the quenched case. Our results are in agreement with 
the indications for a first order phase transition in full QCD with 2 dynamical flavours (in 
the same range of quark masses) obtained in Refs. [11, 16, 17]. 
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Figure 10: F'(/3,Lf,m q ) rescaled by the factor (Lf y. The values of Lf and m q = am are 
displayed in the legend. 



4. Conclusions 

Let us conclude by stressing the main results of this paper. We investigated the nature of 
deconfining phase transition in SU(3) pure gauge theory and in full QCD with two flavors of 
staggered fermions. To locate the phase transition we used the derivative of the monopole 
free energy with respect to the gauge coupling. The monopole free energy is defined by 
means of a gauge invariant thermal partition functional in presence of the abelian monopole 
background field. In the pure gauge case our finite size scaling analysis indicate a weak first 
order phase transition. We get also an estimate of T c /^/a extrapolated to the continuum 
in good agreement with updated value in the literature. Moreover our method has been 
checked in U(l) pure gauge theory giving results in agreement with previous investigations 
which we present in Appendix A. 

In the case of 2 flavors full QCD, we performed simulations by varying spatial lattice 
sizes and quark masses. We find that deconfinement transition in full QCD with 2 degen- 
erate dynamical flavors is consistent with a weak first order phase transition, contrary to 
the expectation of a crossover for not too large quark masses, but in agreement with the 
indications obtained with Refs. [11,16,17]. Our results deserve further investigations that 
we plan to do in the near future. In particular we would like to stress that we have used the 
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Figure 11: U(l) lattice gauge theory. The derivative with respect to (3 of the monopole energy on 
lattices L 4 (L = 16,24,32) and n mon = 20. 



standard pure gauge and staggered fermion action with L t = 4, where finite lattice spacing 
scaling violations could still be important, so that we plan to make use of an improved 
action and/or of a larger value of Lt. We plan to investigate the critical region by means 
of the density spectral method. However, it is worthwhile to stress that, if this result will 
be confirmed the phase diagram of QCD with dynamical flavors should be reconsidered. 
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A. U(l) 

In this Appendix we take into account pure gauge U(l) lattice theory at zero physical 
temperature, essentially for the purpose of showing that using our method the known results 
of compact U(l) are reproduced. We consider U(l) lattice gauge theory in a monopole 
background field. In the continuum the magnetic monopole field with the Dirac string in 




the direction n is 

v ' 2 \x\(\x\ -x-n)' y ' 

where, according to the Dirac quantization condition, n mon is an integer and e is the electric 
charge (magnetic charge = n mon /2e). On the lattice the links belonging to the xt = time- 
slice and to the spatial boundaries of the xt ^ time-slices are constrained as (we choose 
n = x 3 ) 

U$(x) = cos[^ 2 on (*)] +ism[9™\x)] , 



C/| xt (f) = l, 



with 



(X2 ~ X 2 ) 1 



If^(x) = + 



2 | •''mon |^mon| (^3 X;\) (A3) 

(si - Xi) 1 v " ^ 

2 | mon l^-monl (^3 X3) 



In Equation ( |A.3[) (X\, X 2 , X%) are the monopole coordinates and x mon — (x-X). In the 
numerical simulations we put the lattice Dirac monopole at the center of the time slice 
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xt = 0. To avoid the singularity due to the Dirac string we locate the monopole between 
two neighboring sites. 

The abelian monopole condensation is detected by means of £"(/?), the derivative 
of the vacuum energy (with respect to the gauge coupling) in presence of the monopole 



background field and is given by (see Sect. |1.1|) 



E 



1 Z[A exi 
U Z[0] 



(A.4) 



where i?L4 ext ] is the lattice Schrodinger functional with A ext the monopole background 
field and, according to the physical interpretation of the effective action, we denote T in 
Eq. flOD with E. 

Again, to avoid the problem of evaluating partition functions, we compute E'(/3), the 
derivative with respect to (3 of the monopole energy E((3), defined in terms of the lattice 
effective action Eq. Ql.l|) . E'((3) is computed by the analogous of Eq. ( |1.9[ ). 

In Fig. |ll| the data for E'((3), with 

Ji m0 n — 10, on lattices (L — 16, 24, 32) are 
displayed. We observe a sharp peak near the deconfining transition {(3 ~ 1.01). To extract 
quantitative features from the lattice data we do a finite size scaling analysis. To this 



purpose we employ the scaling law Eq. (2.11) to fit the data for E'(f3) on lattices 24 
and 32 4 . The best-fit parameters are reported in Table In Fig. [l^ we report E f (f3, L° s ), 
rescaled by 1/{L C S ) 7 , versus the scaling variable {Lf ) l / u {f3- (3 C ). As we can see the scaling 
law Eq. (|2.11|) holds very well for a wide range of the scaling variable x. By inspecting 



spatial "fixed boundary conditions" (see Section 2.1) 



Ol 



7 



d 



a 



285.2 
±102.0 



2.80 
±0.70 



1.0107 
±0.0010 



0.245 
±0.010 



200.0 
102.8 



1.019 
±0.036 



Table 4: The values of the parameters obtained by fitting Eq. (2.8) to the data for the derivative 



of the monopole energy Eq. (L9) on lattices with volumes 24 , and 32 4 



Table f|we can see that a/v = 4.16±0.22 which is consistent with 7+I = 3.8±0.7. Indeed, 



according to Eq. (A.4), to obtain a sensible result for /3 < [3 C in the thermodynamical limit, 
we must have a/v = 7+1. Indeed, we see that this last relation is satisfied within statistical 
uncertainty. Our results are in good quantitative agreement with previous finite size scaling 
study of the disorder parameter performed in Ref. [36]. Moreover our determination of 
the critical coupling (3 C is consistent with the most recent determination reported in the 
literature [37] . For what concern the order of phase transition our data indicate a weak first 
order phase transition. In fact the parameter v is consistent with 1/d for d = 4. Moreover 
a ~ 1 within errors indicating that the analogous of the disorder parameter of Refs. [8-10] 
(see the discussion at the end of Sect. |2.2|) goes to zero almost continuously as (3 — > C . In 
this sense, we can say that the first order deconfinement transition in U(l) is weaker than 
in SU(3). 
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